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Abstract. - We investigate the electronic structure and properties of Rydberg atoms ex- 
posed to a magnetic quadrupole field. It is shown that the spatial as well as generalized time 
reversal symmetries lead to a two-fold degeneracy of the electronic states in the presence of 
the external field. A delicate interplay between the Coulomb and magnetic interactions in the 
inhomogeneous field leads to an unusual weak field splitting of the energy levels as well as 
complex spatial patterns of the corresponding spin polarization density of individual Rydberg 
states. Remarkably the magnetic quadrupole field induces a permanent electric dipole moment 
of the atom. 



The past two decades have seen substantial progress of our knowledge on highly excited 
Rydberg atoms exposed to homogeneous magnetic fields providing major impact on areas 
such as quantum chaos, semiclassics of nonintegrable systems and properties of magnetized 
structures POEM]. However, so far there exist no investigations on Rydberg atoms in inho- 
mogeneous and/or trapping magnetic field configurations. Apart from being of fundamental 
interest trapped Rydberg atoms have recently been proposed to serve as a tool for quantum 
information processing in mesoscopic atomic ensembles p]. For these applications the atoms 
have to be localized in sufficiently tight traps providing the appropriate confinement of the 
highly excited Rydberg states. Such a tight confinement can be achieved for neutral atoms in 
magnetic traps on atom chips [H] where large field gradients B rs 10 8 ^ are accessible. As a 
prototype example we study here the magnetic quadrupole field [2] that is a key element of 
magnetic trapping We show that Rydberg atoms confined to a quadrupole field possess 
a specific structure and symmetry that lead to two-fold degeneracies of its eigenstates. The 
Rydberg states exhibit unique phenomena such as complex spin polarization patterns and 
magnetic field-induced giant electric dipole moments that can be understood by employing 
the underlying symmetries and analyzing the interplay between the Coulomb and magnetic 
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( 1 )The quadrupole field is not a complete trap in itself 
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interactions. We utilize a one-body approach for the Rydberg atom where the motion of the 
excited outermost electron takes place in the field of a singly positive charged core. The accu- 
racy of this assumption increases with increasing degree of excitation and holds particularly 
well for the frequently used alkali atoms which possess a single valence electron outside a 
closed shell core ( 2 ). We assume that the atomic center of mass (CM) is localized at the cen- 
ter of the quadrupolc field requiring an ultracold CM motion of the atom. The Hamiltonian 
H = Hi + H2 in the presence of the quadrupole field B(r) = B(x, y, —2z) with the vector 
potential A{r) — \ [B{r) x r\ reads 



We have employed spherical coordinates. Ui(i — x, y, z) are the Pauli matrices acting in spin 
space. Here B, Hb are the field gradient and Bohr magneton, respectively. Hi contains the 
Coulomb, paramagnetic and diamagnctic interactions whereas Hi contains the interaction of 
the spin with the magnetic field. Compared to the case of a homogeneous field the Hamilto- 
nian H exhibits a number of major differences. Depending on the value of the field gradient 
B it possesses a strong variability with respect to the appearance of its energy surfaces. The 
paramagnetic term (oc B) is, apart from its proportionality with respect to the angular momen- 
tum L z — j-§^, additionally depending on the z— coordinate and the diamagnetic interaction 
(oc B 2 ) represents an oscillator coupling term of fourth order for the motion perpendicular 
and parallel to the z— axis (see Hi). In addition H2 results in an intricate coupling of the 
spatial and spin electronic degrees of freedom via the inhomogencity of the field. It is possible 
to eliminate the ^—dependence of the Hamiltonian H by applying the unitary transformation 



acting in spin and angular space. This is due to the rotational symmetry associated with the 
conservation of the total angular momentum J z — L z + S z possessing half-integer eigenvalues 
M. This conservation reflects the axial symmetry of the quadrupole field. The eigenfunctions 
of J z (andH) take the appearance I $ m ) = ($ u (r, 9)e^ m - 1 ^, $ d (r, 9)e m ^) with M = m-\,m 
being integer, and where <&"' d are the upper and lower components of the spinor, respectively. 
A close inspection of H yields a further spatial (unitary) symmetry P^,OP z that consists of 
the spatial z— parity operation P z followed by an interchange 0(— cr x ) of the spin components 
and the </>— parity P^ : (f> — > 2ir — <p operation. Additionally H possesses the generalized time 
reversal symmetry TOP z i.e. [TOP z ,H] — 0, where T is the conventional time reversal oper- 
ation. P^OP z and TOP z do not commute with J z but yield e.g. [TOP z , J z ] = {-2J z )TOP z . 
Eigenfunctions to the TOP z -operator are provided by the corresponding linear combination 
|*±) = -1= [|$ m ) ±TOP z \<P m }} with TOP z |*±) = ±|*±). Beyond the above we have the 
additional symmetry operation TP^ that commutes with both the Hamiltonian H and the 




2m e 

H2 = ^bBt (sin 9{a x cos <j) + a y sin <p} — 2o z cos 9) 



+ -^B 2 r A cos 2 9 sin 2 9 



(1) 




(2) 



( 2 )Also, for excited states the spin-orbit and hypcrfine interactions can be neglected due to their rapid drop-off 
with increasing energetical degree of excitation 
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angular momentum J z . Employing {TOP z , J z } — it can be shown that each energy eigen- 
value is doubly degenerate with the two energy eigenstates being |$ m ), |$ m )' = TOP z \& m ) 
i.e. they possess the eigenvalues ±M with respect to J z . The underlying symmetry group is 
Nonabelian and a semi-direct product Coov — Coo (^) C s ■ This is in contrast to the case of a 
homogeneous magnetic field where L z , P z , P,To z P$ constitute the spatial and time reversal 
symmetries, respectively, and form an Abelian group thereby not causing any degeneracies 
due to symmetry. The above symmetries and degeneracies have to be carefully distinguished 
from the two- fold Kramers degeneracy of spin \ systems in the absence of the field. The 
latter is lifted if an external (even homogeneous) field is switched on. The degeneracies found 
here are due to the particular geometry of the quadrupole magnetic field. To investigate the 
electronic structure of the atom exposed to the field in detail we expand the eigenfunctions of 
7i in two-component spinors according to 

*(r,M)= £ w( 4^w-) (M) J ( 3 ) 

where (r, #,</>), YS m ' are spherical coordinates and the spherical harmonics, respectively, and 

Rn' k \r) = J ( n +2k)\ e ~^ T (C r ) k ^n k (C r ) w *th L® being the associated Laguerre polynomials. 
£ is a nonlinear variational parameter to be optimized. We apply the energy variational 
principle linearly optimizing the coefficients c n , - r by solving the corresponding generalized 
eigenvalue problem (GEP). The Hamiltonian and overlap matrices can be calculated analyt- 
ically possessing a band and block structure, respectively, that can be exploited in solving 
the GEP. Our approach to diagonalize the GEP consists of a Krylov-space approach using 
the Arnoldi-decomposition and furthermore applying a shift-and- invert procedure (SJ|U] . This 
allows us to accurately describe highly excited states since the shift-and-invert approach com- 
bined with an optimized value for the parameter ( (1/ y/0.77b\E\ is a good choice) allows 
to converge the eigenvalues in a preselected window of the excitation energy. For the field 
gradients up to B — 10 13 ^ we have studied excitation spectra corresponding to the mag- 
netic quantum numbers M = (— |) — (+|). For the gradients presently accessible for atom 
chips (B < 10 7 — ) all Rydberg states up to n sa 60 (n refers to the corresponding field-free 
principal quantum number which is not a good quantum number in the presence of the field 
but serves as an energetical label) could be calculated with high accuracy. To investigate the 
properties of the excited states in more detail we analyzed the spectrum Ei(B) (i-th energy 
curve), the spatial probability densities W\(r, 0) = r 2 sin#|$ m (r, 9)\ 2 of individual states, the 

z-component of the spatial spin density Ws (V, 0) = | ^[r'elja+l^^'g^j as wen as several 
relevant expectation values such as dipole moments and spin polarizations. Emerging from 
B = we encounter for weak gradients a splitting of the degenerate energy levels that looks 
quite different from what is observed in the case of a homogeneous magnetic field. Focusing 
e.g. onM = i and a given n— manifold the splitting in a homogeneous field takes place into 
two bundles each consisting of almost degenerate n and n — 1 sublevels (see inset of upper 
panel of figure 1), respectively. One of the bundles is (approximately) independent of the field 
and the other one raises linearly with increasing field strength. In the quadrupole field the 
degenerate n— manifold of M = ^—states splits into its 2n — 1 components i.e. energy curves 
Ei(B) (see upper panel of figure 1). Each curve behaves approximately linear with increasing 
field gradient (in the weak gradient regime where the paramagnetic interaction dominates) 
but possesses a different slope. The curves Ei(B) are arranged symmetrically with respect 
to the constant Ei(Q) as can be seen in figure 1. With increasing B the clusters of levels 
widen further and the curves Ei(B) become nonlinear (see middle panel of figure 1). The 
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Fig. 1 - Upper panel: The energy levels as a function of the field gradient B emerging for the 
n = 19 — 23, M — | manifolds. The inter n-manifold mixing for strong gradients is clearly visible. 
Lower right panel: Symmetric linear splitting of the 39 (for E = 0) degenerate energy levels belonging 
to the multiplett n — 20, m = 1 for weak gradients. Inset: The splitting of the same states in 
a homogeneous magnetic field covering the range < B < 0.235 Tesla. Lower right panel: The 
splitting of the same multiplett for a larger range of field gradients covering the regime of Z-mixing. 
Very narrow avoided crossings occur in this regime. 
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Fig. 2 - The spatial probability density Wa (left column, logarithmic representation) together with 
the S'z-density in r, #-space (right column) for the 15-th excited state with M = \ for B = 4.4- lO 9 ^- 
(upper panel) and for the 76— th excited state (lower panel) emerging from n = 9 with M = | for 
the same gradient. The asymmetry of the probability density with respect to the line = 5, i.e. for 
2 — > — z, is clearly visible. 



field gradient for intra n-manifold mixing of different angular momentum i.e. /-states can be 
shown to scale as B cx n~ 6 . In this regime very narrow avoided crossings occur. With further 
increasing field gradient the n-manifolds start to overlap (see lower panel of figure 1) and we 
encounter inter n- manifold mixing that scales according to B cx nT ~ (in a homogeneous field 
the corresponding scaling for inter n-manifold mixing is B cx n~ 2). In this regime the dia- 
magnetic interaction of T~L\ is important and no (not even approximate) symmetries remain. 
Level repulsion and avoided crossings are therefore a characteristic feature of the spectrum. In 
the quadrupole field the spatial probability distributions of the J z — eigenstates are typically 
localized in one half-volume cither above or below the x, y— plane (see left column of figure 2 
which shows the distributions for the 15-th and 76— th excited state for B — 4.4 • 10 9 — ) ( 3 ). 

cm / v / 

With increasing degree of excitation of the eigenstates, the diamagnetic interaction becomes 
important and leads to an additional deformation of the electronic probability density. In a 
homogeneous magnetic field z-parity is a symmetry and Wa is symmetric with respect to re- 
flections at the x, y-plane. Considering the expectation value (r) for the states with increasing 
n the effect of the quadrupole field is twofold: it modifies the distribution of (r) values within 
a single n— manifold and decreases the center of the distribution due to the compression of 
the electronic cloud in the quadrupole field. Looking at the spatial Wg-spin density of the 
electronic states some remarkable properties appear. Figure 2 shows Ws(r,6) for the 15-th 
and 76— th excited state for a field gradient B = 4.4 ■ 10 9 ^. For the 15-th excited state 
Ws shows curved stripes of upward and downward polarized spin, respectively, that match 
with the regions of the localization of the spatial probability density W\. For the 76— th 
excited state a pattern of nested islands appears with each island possessing a certain spin 
polarization and well-localized transition regions separating them. These islands correspond 



( 3 )This property is related to the existence of permanent electric dipole moments as discussed below 
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Fig. 3 - The expectation value of the electric dipole moment along the z— axis for B = 4.4- 10 s -2- for 
the states with M = i. n labels the excitation energy. 



to locally either upward or downward pointing spin and are arranged in a chess board-like 
pattern. The borderlines between the islands correspond to a vanishing ^-component of the 
spin. The intersection of the borderlines i.e. the corners of the islands, represent the nodes of 
the spatial probability densities of the Rydberg states as can be seen from the corresponding 
graph Wa in figure 2. With increasing value of r the shape of the islands become elongated in 
radial direction and finally turn into directed stripes with continuous transitions of the spin 
polarization. The formation of the islands is due to a detailed balance of the interactions 
in the Hamiltonian TL. This does not occur for the atom in a homogeneous magnetic field 
(constant spin polarization) nor for the case of a pure spin in a quadrupole magnetic field 
described by TL2 only. The latter yields a spin polarization density that is independent of r 
showing horizontal stripes in the (r, (^-representation of Ws- These stripes can be found for 
Ws of the Rydberg states in figure 2 for large values of r in a modified form indicating the 
dominance of the spin coupling to the magnetic field in this regime. Analyzing the spin density 
of many states we found that the above behaviour is generic. The spin polarization patterns 
are features of individual electronic eigenstates and become increasingly more detailed with 
increasing degree of excitation of the state considered. Transitions among eigenstates to the 
total angular momentum J z obey the selection rules AM = and A M = ±1 for dipole 
transitions via linear and circular polarized light, as it is the case without the presence of the 
field. Linearly z— polarized transitions between TOP z -eigenstates must involve states with 
different TOP z symmetries in order to possess a nonvanishing dipole strength. The fact that 
the quadrupole field causes a nonsymmetric charge distribution with respect to the horizontal 
x, y— plane leads to the following peculiar properties of the atom: Electronic eigenstates to J z 
possess a nonvanishing permanent electric dipole moment e(z) only along the symmetry axis 
of the quadrupole field i.e. the external magnetic field induces a permanent electric dipole 
moment of the atom. Figure 3 illustrates the distribution of dipole moments for a variety of 
states belonging to the manifolds n < 65. The variance of the distribution of dipole moments 
increases strongly with increasing degree of excitation n thereby showing a transition from 
a regular alignment to an irregular spreading. In a homogeneous magnetic field the defor- 
mation of the charge distribution is (due to parity symmetry) such that the electric dipole 
moment vanishes. One possibility to probe the above-described properties of Rydberg atoms 
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in quadrupole magnetic fields would be to perform spectroscopy of single atoms in traps on, 
preferably, atom chips since these possess currently the strongest available field gradients. 
This would provide us with detailed information on the level splitting and evolution with 
increasing degree of excitation. 

Already in the presence of a homogeneous magnetic field it is well-known that the CM and 
electronic motion of an atom do not decouple |1UU11|IT2] . To enter the corresponding regime 
where the residual coupling becomes important certain parameter values (excitation energy, 
CM energy etc.) have to be addressed. A variety of intriguing phenomena due to the mixing 
of the internal and CM motion such as the diffusion of the CM or giant dipole states are then 
observed |13U14lH5| . In the quadrupole field the assumption that the atom is ultracold will 
certainly minimize the CM motional effects. Nevertheless, a residual coupling is unavoidable 
and its impact on the electronic structure is, at this point, simply unkown: A full treatment of 
the two-body system certainly goes beyond the scope of the present investigation and requires 
both from the conceptual as well as computational point of view major investigations. On the 
other hand one should note that the symmetries discussed here equally hold for the moving 
two-body system i.e. the total angular momentum is conserved and the unitary as well as 
antiunitary spin-spatial symmetries, now applied to both particles, are also present. 

Beyond the above, one can speculate about potential applications of the magnetic field- 
induced permanent electric dipole moments of the atoms. Populating with a laser excited 
states with a desired dipole moment for certain atoms within an array of single atom traps 
can open the route to a controlled interaction between the atoms which is currently of major 
interest for quantum information processing |16U17llolll8| . Discussions with Ofir Alon and M. 
Anderson are gratefully acknowledged. J.S. acknowledges financial support by the European 
Union contract numbers IST-2001-38863 (ACQP). 
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